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Salmhofer [Commun. Math. Phys. 194, 249 (1998)] has recently developed a new renormaliza- 
tion group method for interacting Fermi systems, where the complete flow from the bare action of 
a microscopic model to the effective low-energy action, as a function of a continuously decreasing 
infrared cutoff, is given by a differential flow equation which is local in the flow parameter. We 
apply this approach to the repulsive two-dimensional Hubbard model with nearest and next-nearest 
neighbor hopping amplitudes. The flow equation for the effective interaction is evaluated numeri- 
cally on 1-loop level. The effective interactions diverge at a flnite energy scale which is exponentially 
small for small bare interactions. To analyze the nature of the instabilities signalled by the diverg- 
ing interactions we extend Salmhofers renormalization group for the calculation of susceptibilities. 
We compute the singlet superconducting susceptibilities for various pairing symmetries and also 
charge and spin density susceptibilities. Depending on the choice of the model parameters (hop- 
ping amplitudes, interaction strength and band-filling) we find commensurate and incommensurate 
antiferromagnetic instabilities or d-wave superconductivity as leading instability. We present the 
resulting phase diagram in the vicinity of half-fllling and also results for the density dependence of 
the critical energy scale. 
PACS: 71.10.Fd, 71.10.-w, 74.20.Mn 
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I. INTRODUCTION 

One of the striking aspects of the high-temperature 
superconducting cuprates is the sensitive dependence of 
their physical properties on the charge carrier concen- 
tration in the copper-oxide planes, which can be con- 
tinuously varied by doping the inter-plane region. In 
the doping-temperature phase-diagram one generically 
finds an antiferromagnetic insulator and a superconduct- 
ing phase with d-wave symmetry, with strongly doping 
dependent transition temperatures in each case . 

The two-dimensional Hubbard model Q| is a promising 
prototype model for the electronic degrees of freedom in 
the copper-oxide planes. It has an antiferromagnetically 
ordered ground state at half-filling and is expected to be- 
come a d-wave superconductor at moderate doping away 
from half-filling 

Although the Coulomb interaction in the cuprates is 
certainly rather strong, there has been considerable re- 
cent interest in the weak coupling sector of the 2D Hub- 
bard model. Besides the applicability of (semi-) analyti- 
cal methods at weak coupling and the general experience 
that many strongly interacting systems are more or less 
continuously connected to corresponding weak coupling 
systems, a major reason for this interest is that even the 
weakly interacting 2D Hubbard model exhibits an ex- 
traordinarily rich behavior as a function of the carrier 
density and other model parameters. Conventional per- 
turbation theory breaks down for densities close to half- 
filling, where numerous competing infrared divergences 
appear as a consequence of Fermi surface nesting and 
van Hove singularities. These divergencies can in prin- 
ciple be treated by suitable self-consistent resummations 
of perturbative contributions to all orders in the cou- 



pling constant. Most notably the so-called fluctuation- 
exchange approximation Q turned out to be able to de- 
scribe various expected physical properties. However, 
a completely unbiased selection of Feynman diagrams 
that takes into account all possible particle-particle and 
particle-hole channels on equal footing would require at 
least the self-consistent summation of all parquet dia- 
grams which is still beyond present numerical pos- 
sibilities for sufficiently large systems and low tempera- 
tures. 

Renormalization group (RG) methods are presently 
the most promising and best controlled approach to low 
dimensional Fermi systems with competing singularities 
at weak coupling. Such methods have been developed 
long ago for one-dimensional systems where, combined 
with exact solutions of fixed point models, they have been 
a major source of physical insight Early RG studies 
of two-dimensional systems started with simple but inge- 
nious scaling approaches to the 2D Hubbard model, very 
shortly after the discovery of high- Tc superconductivity 
|p|-|Tot. These studies focussed on dominant scattering 
processes between van Hove points in k-space, for which 
a small number of running couplings could be defined and 
computed on 1-loop level. Spin-density and supercon- 
ducting instabilities where identified from divergencies 
of the corresponding correlation functions. Recently, the 
early scaling approaches have been revisited by various 
authors to extract further physical properties, such as a 
possible pinning of the Fermi level at the van Hove singu- 
larity ||ll[] , extended saddle points 1 11 ij] , and a possible 
gap formation on parts of the Fermi surface near the van 
Hove points [ p^ . Scaling theories with few running cou- 
plings have also been used to analyze instabilities associ- 
ated with flat Fermi surface pieces ||l^,|l5| and inflection 
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points on the Fermi surface . 

A major complication in two-dimensional systems 
compared to ID is that the effective interactions cannot 
be parameterized accurately by a small number of run- 
ning couplings, even if irrelevant momentum and energy 
dependences are neglected, since the tangential momen- 
tum dependence of effective interactions along the Fermi 
"surface" (a line in 2D) is strong and important in the 
low-energy limit. This has been demonstrated in partic- 
ular in a 1-loop RG study for a model system with two 
parallel flat Fermi surface pieces In an impressive 

series of recent papers Zanchi and Schulz have de- 
veloped a new renormalization group approach for inter- 
acting Fermi systems, which is based on a flow equation 
derived by Polchinski in the context of local quan- 
tum field theory. In this RG version the complete fiow 
from the bare action of an arbitrary microscopic model to 
the effective low-energy action, as a function of a contin- 
uously decreasing infrared energy cutoff, is given by an 
exact differential flow equation. Zanchi and Schulz have 
applied this approach to the 2D Hubbard model (with 
nearest neighbor hopping) in a 1-loop approximation, 
with a discretized tangential momentum dependence of 
the effective 2-particle interaction. The presence of anti- 
ferromagnetism and d-wave superconductivity as major 
instabilities of the model close to half-filling was thereby 
confirmed. 

The development of continuous renormalization group 
methods for interacting Fermi systems has made further 
progress with a most recent work by Salmhofer [po[ . By 
expanding the effective action in Wick-ordered monomi- 
als instead of bare monomials he obtained an exact fiow 
equation for the effective interactions with a particularly 
convenient structure: The /3-function is bilinear in the ef- 
fective interactions and local in the flow parameter, i.e. it 
does not depend on effective interactions at higher energy 
scales. 

In this work we apply Salmhofer's RG version to the 
two-dimensional Hubbard model with nearest and also 
next-nearest neighbor hopping amplitudes, concentrat- 
ing on the most interesting electron density regime near 
half-filling. We evaluate the flow of effective (2-particle) 
interactions on 1-loop level, neglecting the irrelevant en- 
ergy dependence and also the irrelevant normal momen- 
tum dependence, but keeping the important tangential 
momentum dependence. As in previous RG calculations, 
the effective interactions diverge at a finite energy scale, 
which is exponentially small for a small bare interac- 
tion. To analyze the physical nature of the instabil- 
ities signalled by the diverging interactions we extend 
Salmhofer's RG version for the calculation of susceptibili- 
ties. We compute charge and spin susceptibilities and sin- 
glet superconducting susceptibilities for various pairing 
symmetries. Depending on the choice of the model pa- 
rameters, hopping amplitudes, interaction strength and 
band-filling, we find commensurate or incommensurate 
antiferromagnetic instabilities or d-wave superconductiv- 
ity as leading instability, in qualitative agreement with 



previous calculations with other RG versions. We present 
the resulting phase diagram of the two-dimensional Hub- 
bard model near half-filling and results showing how the 
critical energy scale decreases away from the half-filled 
perfect nesting case. 

We finally note that powerful renormalization group 
techniques with a discrete successive reduction of the in- 
frared cutoff have recently opened the way towards a rig- 
orous non-perturbative control of interacting Fermi sys- 
tems for sufficiently small yet finite coupling strength 
Significant rigorous results have already been de- 
rived for 2D systems These mathematical works 
show in particular that all weak coupling instabilities in 
interacting Fermi systems can be obtained systematically 
from a renormalization group analysis. 

This article is organized as follows. In Section II 
we briefly review Salmhofer's RG for general interact- 
ing Fermi systems, present the explicit fiow equations 
for effective 2-particle interactions on 1-loop level, and 
also derive the flow equations for susceptibilities. In Sec- 
tion III we first discuss how the 1-loop flow equations are 
solved, and then present numerical results for the Hub- 
bard model with a discussion of their physical content. 
Section IV closes the presentation with a short summary 
of major results and some ideas for further developments. 

II. RENORMALIZATION GROUP EQUATIONS 

In this section we review Salmhofer's renormalization 
group approach for general interacting Fermi systems, 
present the explicit flow equations for effective 2-particle 
interactions on 1-loop level, and finally derive 1-loop fiow 
equations for several susceptibilities, which will later be 
used for our stability analysis of the 2D Hubbard model. 

A. Functional integral representation 

We consider a system of interacting spin-i fermions 
with a single particle basis given by states with a (crys- 
tal) momentum k, a spin projection a G {Tii}5 ^-nd ^ 
kinetic energy Ck- The properties of the system are de- 
termined by an action 

K 

where K = (fco,k, ct) is a multi-index containing the 
Matsubara frequency fco in addition to the single par- 
ticle quantum numbers; ipK and tpK are Grassmann 
variables associated with creation and annihilation op- 
erators, = — is the single-particle energy rel- 
ative to the chemical potential, and Vl^p^ip] is an arbi- 
trary many-body interaction. The non-interacting single- 
particle propagator of the system is given by 

^(^) = (2) 

iko - 4k 



2 



All connected Green functions can be obtained from the 
generating functional 

g[rj,7]] = logj / d/.c[V'>]e-^['^''^leW'-'')+('^''^4 (3) 



with the normalized Gaussian measure 

dfic[ip,i'] = 



,('/',c-V) 



K 



K 



(4) 



Here we have introduced the short-hand notation 
(X,^) = EkXk^k and iC-'^)K = C-\K)^k for 
arbitrary Grassmann variables xk and ipK- Note the 
identity 



which implies that g[ri,fj] = —{fj,Cri) in the non- 
interacting case Vltp, ip] — 0. The connected m-particle 
Green functions are given by 



Gm{K[, . . . ,K'j-^; Ki, . . . , K„i) — 



dr]K[ ■ ■ ■ driK^ dr]K^ . . . drjK^ 



ll-'Ki)c 



r}—ri—0 

(6) 



where {■ ■ ■)c is the connected average of the product of 
Grassmann variables between the brackets. 

The renormalization group equations are most conve- 
niently formulated for another generating functional, the 
effective interaction |24 



v[x, x] = - log / d^lc[i^, ^] e-^IV'+x.^+x] 



(7) 



A simple substitution of variables yields the relation 

-V[x, X] = iv, Cv) + G[V, V] where x = Cv, Cv 

(8) 



Hence, functional derivatives of V[x,x] gen- 
erate connected Green functions divided by 
C{Ki) . . . C{Kra) C{K[) . . . C{KIJ, i.e. (non-interacting) 
propagators are amputated from external legs in the cor- 
responding Feynman diagrams. The term [fj, Crj) cancels 
the non-interacting part of Q[f]^i]\ such that V[x,x] = 
for V'] = 0- Hence, the non-interacting propaga- 
tor is subtracted from the one-particle Green function 
generated by V. 



B. Wick-ordered continuous RG 

We now briefly review the derivation of a continuous 
renormalization group equation for Wick-ordered (ampu- 
tated) Green functions, as first derived in the context of 
interacting Fermi systems by Salmhofer [ pO[ |. 

The original system is endowed with an infrared cutoff 
at an energy scale A > by replacing the bare propagator 
C{K) with 



ikQ - 



(9) 



Here e^{K) is a function that vanishes for y/k^ + ^ 

A and tends to one for •\//co + Ck ^ ^- ^^^^ ^^"^ i^" 
frared singularity of the propagator at fco = and = 
(corresponding to the non-interacting Fermi surface in k- 
space) is cut off at scale A. A simple choice for <dy{K), 
which we will later adopt in our numerical calculations. 



©dCkl-A) 



(10) 



where Q is the usual step function. With this choice 
single-particle states close to the Fermi surface are 
strictly excluded from the theory. Alternatively, one may 
also use a smooth cutoff function. 

The generating functional for connected Green func- 
tions and the effective interaction constructed with 
(instead of C) will be denoted by g^[j],fj] and V"^[x,x]: 
respectively. The original functionals Q and V are recov- 
ered in the limit A — > 0. It is not hard to show that the 
effective interaction satisfies the following exact renor- 
malization group equation |20|,|25|| 



where C^{K) = dC^{K)/dK. With the initial condition 
V^nx,x] = ^[X,X] (12) 

this equation determines the flow of uniquely for all 
A < Aq. The initial value Aq must be chosen such that 
Qy{K) = for all K and A > Aq. For the step function 
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cutoff introduced above, Aq is the maximal value of |^k|- 

An expansion of the functional V^[x,x] in the renor- 
malization group equation ([ll|) in powers of xk and xk 
would lead to Polchinski's [ |19[ flow equations for am- 
putated connected Green functions, which have been ap- 
plied to the 2D Hubbard model by Zanchi and Schulz . 
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Alternatively, one can also expand with respect to Wick- 
ordered monomials 



1 



m=0 ^ ' K 



E E 

.....ifm K[,...,K'^^ 



xV^iK[,...,K'^;Ki,...,K„,)e 



-A, 



J = l 

(13) 



The exponent in the Wick-ordering operator is the dif- 
ferential operator 



92 



dxK dxK 
with the propagator 

D^{K) ^ C{K)-C^{K) 



iko - 



(14) 



(15) 



where &'^{K) = 1-Q^{K). Note that Z?^ contributes in 
the infrared regime excluded from . The Wick-ordered 
monomials reduce to pure monomials in the limit A — » 0, 
since D^{K) ^ in that limit. Hence, the functions 
tend to the usual amputated connected Green functions 
for A ^ 0. 



Inserting the expansion (113) 
hand side of the RG equation 

the A-derivation acting on the coefficients Vj^ or on the 
Wick-ordered monomials, respectively. Since 



for V"^[x, x] on the left 
11]) yields two terms, with 



d_ 
OA 



-A„A 



YIXK'^XK, = 

K 



dXK dxK 



WXK'^XK, 

(16) 



and t)^ = — C"^, the second term on the left hand side 
cancels the first term on the right hand side of ( pT| ) . Only 
the second term, quadratic in V^, remains. Expanding 
this term with respect to Wick monomials and compar- 
ing coefficients, one can express the A-derivative of 
as a bilinear form of all the other functions ||2^. A 
graphical representation of these flow equations is shown 
in Fig. 1. 

The precise general equation for has rather compli- 
cated combinatorial factors and will not be written here, 
since we will compute only V.^ in a 1-loop approxima- 
tion. Note that one internal line in Fig. 1 corresponds 
to and the others to the low energy propagator D^. 
Note also that the flow equations are local in A. 



C. 1-loop equations 

To detect dominant instabilities of the system in the 
weak coupling limit, it is sufficient to truncate the infinite 
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FIG. 1. Diagrammatic representation of the flow equation 
for the Wick-ordered amputated Green-functions Vj^. The 
internal line with a slash corresponds to t)^ , the others to 
D^; all possible pairings leaving m ingoing and m outgoing 
external legs have to be summed. 
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FIG. 2. Flow equation for the vertex function F^ in 1-loop 
approximation with the particle-particle channel (BCS) and 
the two particle-hole channels (ZS and ZS'). 



hierarchy of flow equations described by Fig. 1 at one- 
loop level, and neglect all components of the effective 
interaction except the two-particle interaction V.^. Since 
self-energy corrections are also neglected, reduces to 
the one-particle irreducible two-particle vertex Junction^ 
and will therefore be denoted by F^ from now on. Sum- 
ming all possible pairings of two vertices in a one-loop 
diagram as in Fig. 2, one obtains the flow equation 



^Elr[^^^)^^(^')] 



I3V ^ dK 

^ K,K' 

- T^{K[,K; K,, K') F^(if', K'^-K, K^) 

+ T^iK'^, K- Ki, K') T^{K\ K[;K, K2) 

(17) 



where [3 is the inverse temperature and V is the volume 
of the system. The three terms on the right hand side 
are the contributions from the Cooper (or BCS) chan- 
nel and the two zero-sound channels (ZS and ZS'). Note 
that for translation invariant systems momentum con- 
servation implies that T^{K[^ Ki^] Ki, K2) ^ only if 
ki -\- k2 = k'l + k'2, so that the sum over k and k' in ( p^ ) 
is reduced to a single energy-momentum sum. 



4 



For a spin-rotation invariant system the spin structure 
of the vertex function can be written as 



~^ (^1' ^2' ^1' ^2) TcrJjO-j;. 



0-2 I 

cri,(T2 



where 



(7, ,(7n',<Tl ,(72 



(18) 



(19) 



are the projection operators on singlet and triplet states 
in a two-particle spin space, respectively. The antisym- 
metry of with respect to Ki -s-^ K2 or K{ K2 
implies that is symmetric and antisymmetric un- 
der exchange of the variables ki and ^2 or k[ and 
Carrying out the spin sum in the flow equation one ob- 
tains 

5AF^(fc;,fc^;fci,fc2) = 

E E ra'/3r'(^i'fc2;fci,fc2) + 

i—s,t j—s,t 

+ Cl%f3f{k{,k'2\ kiM) + Cl%pf{k[,k'2: h, k2)] 

(20) 

for a = s,t, where the coefficients Caij can be grouped 
in matrices 



r-BCS _ fi o\ ^zs _ r'zs' _ 1 / 

- (0 0) - - 4 I 



-1 3^ 
3 3y 



CBCS = (0 0^ C7f =Cf' = i(}D 
and the "/3-functions" are given by 
(3^^^k[,k'2;k^,k2)= ^^aA[i?^(fc)D^(fc')] X 



(21) 



fc,fc' 



xr^ik[,k'2;k,k')Tf{k,k';kuk2) 
0ff{k[,k'2-MM) = -;^E^A [i?^(fc)D^(A:')] X 



fe.fe' 



xrf(fc'i,fc;fci,fc')r^'^(fc',fc^;fc,fc2) 
I3ff {k[,k'2-M,k2) = -f3!f{k'2,k[;ki,k2) 

(22) 

We finally list some useful relations for the vertex func- 
tion following from general symmetries. Time-reversal 
invariance implies 

r''{K[,K'2;K,,K2) =a[a'2aia2 x 

X T^{nKi,nK2; nK[,nK^) 

(23) 

where TZK = (A'o,— k, — cr) for K = (fco;k, cr), and the 
number 1 (—1) is assigned to a = f (|) in the pref ac- 
tor. Assuming in addition spatial reflection invariance 
and spin rotation invariance, one obtains 

r^{K[, K'2, K,,K2) = r^(i^i, K2; K[, K'2) (24) 



Prom the behavior under complex conjugation 

f ^(i^;, K'2, Ki,K2) = r^(^i, K2; K[, K'2) (25) 

with K = (— /cq, k, cr) one can then deduce that the ver- 
tex function F^(i^(, ii'2; Ki,K2) is real, if all the energy 
variables vanish. 

D. Susceptibilities 

To identify possible instabilities of the system we com- 
pute various susceptibilities, i.e. the linear response of 
the system to small external fields. 

Application of an external field h leads to an additional 
contribution to the action 

V'[h; V>, Vi] = - E W)0{q) + h{q)0{q)\ (26) 

where 0{q) is a bilinear form in the Grassmann vari- 
ables, 0{q) its hermitian conjugate, and h{q) is the com- 
plex conjugate of h{q). We will compute the response to 
fields coupling to the charge density 



(27) 



and the spin density in ^;-direction 

s'^iq) = E [^fc-9,TV'fe,T - V'fe-g,iV'fc,i] (28) 



and also the response to pairing fields coupling to the 
singlet pair operator 



A{q) = ^ rf(k + q/2) Vfe+«,T^-fc, 



(29) 



where <i(k) is a function with even parity specifying the 
orbital symmetry of the pairing operator (s-wave, d-wave 
etc.). The linear response of the expectation value 



/ \ w ^ d\ogZ[h] 

m{q) = V {Oiq)) = y-^ 



with the partition function 

z[h] = Jdi^c[i',i'] 

is given by the susceptibility 
dm{q) 



-V[^,ip]-V'[h;ip,ip] 



x{q) = 



dh{q) 



= V-^0{q)0{q))l 



h=0 



h=0 



(30) 



(31) 



(32) 



We consider only systems which are translation invari- 
ant, spin-rotation invariant, and charge conserving in 
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the absence of the external field h. In the normal (not 
symmetry-broken) phase the expectation values {0{q)) 
then vanish for h{q) 0, except the expectation value 
{p{0)), which yields the average particle number of the 
system. 

The renormalization group equation ( pj]) is not af- 
fected by the presence of the additional term V' in the 
bare action, since an arbitrary many-body interaction 
was allowed anyway. Only the initial condition of the 



flow is modified to 

V'"'[h;x.x]^V[x,x]+V'[h;x,x] (33) 

which leads to an ^.-dependent effective interaction 
V^[h; X, x]- When a pairing field is coupled to the system, 
the expansion of V^[h; x, x] with respect to Wick-ordered 
monomials will also contain monomials where the num- 
ber of creation and annihilation variables is not equal, 
and Eq. (|l^) must be generalized to 



V%;x,x]= E 



7 't7 ' 



mm. 

7n,Ti— K 



E 



E 

K' Ki,...,K„ 



V±{[h];K[, ...,K'„^;K,,..., K,,) e'^o^ xk; ■ ■ ■ Xk'^ Xk^ ■ ■ ■ Xk, (34) 



To obtain the linear response of the system, we expand 
the effective m-body interactions Vj^ (or its generaliza- 
tion 1^„) in powers of h. The effective 0-body interaction 
can be expanded as 



Vn^h] = V^^ ^vY^x^iq) m Kq) + 0{h^ 



(35) 



Note that V^[h] converges to the grand-canonical poten- 
tial = — log Z [ft,] for A 0, and x^il) converges to 
the susceptibility x{q)- If ^ couples to the charge density, 
there is also a linear contribution —N^[h{Q) + h{Q)\ in the 
expansion of V,^ [h\ , where converges to the average 
particle number for A ^ 0. For a field coupling to charge 
or spin density, the effective 1-body interaction becomes 

V,^{[h];K';K)^V,^{K';K)- 

-Y,[R^{q-K'-K)h{q) + h.C.]+0(ft2) 

(36) 

The spin structure of the renormalized vertex is 

R^{q-K'-K)^5„,,RUq-:k'-k) (37) 
for the charge vertex and 



R^{q;K';K)=a6^,„R^{q;k';k) 



(38) 



for the spin vertex. For a pairing field, V-/^ has only 
quadratic terms in h, but the off-diagonal effective inter- 
actions 

Fo^2(H; i^i, if2) = -2 E ^^(9; Ki,K2) h{q) + 0{h') 
V^%i[h];K[,K^) = -2 E ^""(9; K^K!,) h{q) + 0{h^) 

(39) 

have linear (and higher order) contributions. The spin 
structure of the renormalized singlet pairing vertex is 

i?^(q;i^i,i^2) =CTi5.i,-..i?^(g;fci,fc2) (40) 



Effective two- and many-body interactions can be ex- 
panded similarly. Inserting the expansions of the effective 
interactions into the flow equations and comparing 
the coefficients of contributions with equal powers in the 
external field, one obtains fiow equations for these coef- 
ficients, especially for x^ £ind R^. 

We will again truncate the infinite hierarchy of flow 
equations at 1-loop level, and neglect (zero-field) self- 
energy terms by setting V^{K']K) = 0. To obtain the 
linear response susceptibility it is then sufficient to solve 
the two flow equations for the susceptibility x^ ^'Hd the 
renormalized charge, spin, or pairing vertex R^ repre- 
sented diagrammatically in Fig. 3. If ft couples to the 
charge or spin density, these equations read 



Q^X v^y 



('^) =7^ E dK[DHK)D\K')\ X 

(41) 

X R^{q- K'- K) R^i-q; K; K') 



and 



±^R\q;K[;K,) ^ - ^ ^ 9a [DHK)D\K')\ 



K,K' 



X R^{q-K'-K)T''{K[,K-Ki,K') 



(42) 



with the initial condition R^° {q; K' ; K) = 5a' 
for the charge and R^°{q;K';K) — crSa'aSk-k' ,q 



k — k' ,g 

for the 

spin vertex part. Note that R^{q; K'; K) ^ only if 
k' = k — q. If ft couples to pair creation and annihilation 
operators, the flow equations become 



§X^i<D=-^ E d^[D\K,)D\K,)] X 



kTk. ^ ' (43) 

X R^{q;Ki,K2)R^{q;Ki,K2) 
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FIG. 3. Flow equations for (a) the susceptibilities 
and (b) the response vertices in 1-loop approximation. 

and 

(44) 

with the initial condition {q; Ki, K2) — 

cri5ai-<j'Jki+k2.qd{ ^^2^^ ). Note that 
R^{q]Ki,K2) ^ only if fci + fc2 = q- The initial 
condition for the susceptibility is x^"{<l) = in all cases. 

The vertex function is obtained by solving the flow 
equations in the absence of a perturbing field. One 
can then solve the linear differential equation for R^ and 
finally integrate the flow equation for x^- Note that for 
the special case of a flat fermi surface these 1-loop equa- 
tions have the same structure as the so-called fast parquet 
equations but here the flow variable is the cutoff A 
instead of an external energy or momentum variable. 

III. APPLICATION TO THE 2D HUBBARD 
MODEL 

We will now apply the general renormalization group 
method derived in the previous section to the two- 
dimensional Hubbard model, the main aim being an anal- 
ysis of the leading instabilities of the system at weak cou- 
pling. 

A. Hubbard model and Fermi surface 

The Hubbard model is a lattice electron model with 
the simple Hamiltonian 

H = J2Y.^ijclci, + uY,rii^nj^ , (45) 

ij j 

where and Cio- are the usual creation and annihilation 
operators for fermions with spin projection a G {f, |} on 
a lattice site i. We consider the Hubbard model with a re- 
pulsive interaction L/ > on a (two-dimensional) square 
lattice with a hopping amplitude 

{—t if i and j are nearest neighbors 
— t' if i and j are next-nearest neighbors (46) 
else 




FIG. 4. The Fermi surfaces of the non-interacting 2D Hub- 
bard model with t' = (a) and t' = — 0.16Z; (b) for various 
choices of the chemical potential jj,. 

The corresponding dispersion relation of non-interacting 
particles is 



Ek — — 2i(cos kx + cos ky) — 4t'(cos cos ky) (47) 



This dispersion relation has saddle points at k = (0, tt) 
and (tt, O), which lead to logarithmic van Hove singulari- 
ties in the non-interacting density of states at the energy 
e = ivH = 4^'. In Fig. 4 we show the Fermi surfaces 
of the non-interacting system for various electron densi- 
ties n for the case without next- nearest neighbor hopping 
(t' — 0) and a case with t' < 0. 

For a chemical potential = evH the Fermi surface 
contains the van Hove points. In this case a perturba- 
tive calculation of the two-particle vertex function leads 
to several infrared divergencies already at second order 
in U, i.e. at 1-loop level (see the Feynman diagrams in 
Fig. 2) P~p^. In particular, the particle-particle channel 
diverges as log^ for vanishing total momentum ki + k2 , 
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and logarithmically for ki + k2 = (tt, tt). The particle- 
hole channel diverges logarithmically for vanishing mo- 
mentum transfer; for momentum transfer (tt, tt) it di- 
verges logarithmically if <' 7^ and as log^ in the special 
case t' — 0. Note that /i = for t' — corresponds to 
half-filling {n = 1). For t' = there are also logarith- 
mic divergences for all momentum transfers parallel to 
(tt, tt) or (tt, — tt) due to the strong nesting of the square 
shaped Fermi surface. For ^ evH only the usual log- 
arithmic Cooper singularity at zero total momentum in 
the particle-particle channel remains. However, the ad- 
ditional singularities at /i = EvH fead clearly to largely 
enhanced contributions for small /i — Evh, especially if t' 
is also small. 

Hence, for small /i — evH one has to deal with compet- 
ing divergencies in different channels. This problem can 
be treated systematically by the RG method described 
in Sec. II. 



B. Parameterization of the vertex function 

We now prepare for a computation of the 1-loop flow of 
the vertex function and susceptibilities. The flow equa- 
tions derived in Sec. II cannot be solved analytically. 
Even a numerical solution is not possible with the full 
energy and momentum dependence of the vertex func- 
tion, which must therefore be suitably simplified. 

The dependence of the vertex function on energy vari- 
ables can be neglected completely without much damage, 
because it is absent in the bare interaction, and irrelevant 
(in the sense of power counting) in the low-energy limit 
(see, for example, Ref. |26|). We therefore approximate 



(48) 



Choosing an energy independent cutoff function 8^ (k) , 
the Matsubara sums on the right hand side of the flow 
equations can then be performed analytically. One thus 
obtains 



5Ar^(k'i,k'2;ki,k2 



EE 

i—s,t j—s,t 



/3B.cS(k;, k^; ki, k2) + C^f,/3ff (kl, k^; ki, k2) 



(49) 



for a = s, i, where the /3-functions are now frequency independent and read 

S(k;,k^;ki,k2) =^E^A [e^(k)e^(k')] ^^^!^i--^rf(k;,k'2;k,k')r,^(k,k';ki,k2) 
A^S(k;,k'2;ki,k2) = - i^9A [e^(k)e^(k')] ^^%^^|^rf(k;,k;ki,k')r,^(k',k^;k,k2) (^o) 

Pf (k'l , k'2 ; ki , k2 ) = - /3f/ (k'2 , k'l ; ki , k2 ) 



with the Fermi function /(^) — \e^^ + l] . Note that 
momentum conservation implies that k and k' are re- 
lated by k + k' = ki + k2 in the Cooper channel and by 
k-l-k'^ = k'-f ki in the zero sound channel, such that only 
one independent momentum variable has to be summed. 
For a step cutoff function 0:^(k) = 8(A — |^k|) one has 

[e^(k) (k')] = <5(A - i^kl) e(A - |^k' I) + k ^ k' 

(51) 

such that the k-sum can be reduced to a one-dimensional 
integral in the thermodynamic limit. To this end we sub- 
stitute kx and ky by the new variables ^ = ^k and the 
angle (j) between k and the x-axis in k-space, and use 

(52) 



I 

where J{^,(f)) is the Jacobian associated with the trans- 
formation of variables. Since the integrand contains a 
factor (5(A — |^|), the ^-integration can be performed by 
hand, leaving only the angular integral over to be done 
numerically. 

The flow equation can be solved only if also the mo- 
mentum dependence of the vertex function is simplified. 
At least for weak coupling (in practice also for moderate 
ones), the vertex function acquires strong momentum de- 
pendences only for momenta close to the Fermi surface. 
Note that for the Hubbard model the bare vertex func- 
tion F"^" does not depend on momentum at all. Instabil- 
ities are signalled by divergencies of the vertex function 
F^ for momenta close to the Fermi surface and small A. 
Hence we will focus on the flow of the vertex function 
with momenta close to the Fermi surface. 

The intermediate momenta k and k' are constrained 
to a momentum shell around the Fermi surface, since 
ICkUCk'l < A (see Eq. (^). Hence, the values of the 
vertex function at momenta within that A-shell govern 
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their own flow. 

For fixed finite A, the dependence of r^(k']^, kj; ki, k2) 
on ^ki etc. is regularized by the cutoff for |^ki |, • • ■ < A. 
For momenta within the A-shell one may therefore ap- 
proximate the vertex function by 

r^(ki,k2;ki,k2) « r^(k^i,k_Fi + k_F2 - k'^^; k_Fi, k_F2) 

(53) 

where k^^i etc. are projections of ki etc. on the Fermi 
surface (see Fig. 5). Note that strong momentum de- 
pendences of the effective vertex are built up only by 
contributions with intermediate momenta k and k' (on 
the right hand side of the flow equations) which are close 
to the Fermi surface, because for such momenta the frac- 
tions f^^f^ll^J^^k') ^g^yj i^e big. Hence, for the most im- 
portant momenta, the error made by the projection is 
relatively small (even if A is not small), because these 
momenta are close to their projected counterparts. The 
projected vertex function can be parameterized by three 
angles 4>i,4'2,4'3 associated with ki^i, ki?2 and k^j^, re- 
spectively, i.e. 

(54) 

The angular dependence turns out to be strong for small 
A and cannot be neglected. 

With the above projection procedure only functional 
dependences which are irrelevant in the low-energy limit 
have been neglected (see, for example, Ref. |^). The 
approximation (^3|) is asymptotically exact for A — > 
(for momenta within the A-shell) and, for the Hubbard 
model, also for A = Aq. At intermediate stages of the 
flow there are of course corrections. To assess their im- 
portance, we have refined the parameterization of the 
vertex function in some test runs by using a second sur- 
face of constant energy in momentum space as a target 
for momentum projections. In addition to the Fermi sur- 
face a canonical choice is the van Hove surface defined 
by Ck = EvH, where EvH is the energy at which the den- 
sity of states has a van Hove singularity due to saddle 
points of Ek. Momenta are then projected on that sur- 
face (Fermi or van Hove) which is energetically closer. In 
this way scattering processes between van Hove points, 
which are particularly important at scales A — |evH — mIi 
are treated more accurately. In addition to the three an- 
gles 4>i,4>2, (t>3 one needs three binary variables 1^1,^2, 
to label the closest projection surface for ki, k2 and k'^. 

The parameterization of the vertex parts required 
for the computation of susceptibilities is done in a sim- 
ilar fashion. We concentrate on static susceptibilities 
(go = 0) and neglect the (irrelevant) energy dependence 
of . For the charge and spin density vertex we approx- 
imate 

i?^,s(q;k';k) «i?^^s(q;k;^-q;kj.) (55) 

where kp is the projection of k on the Fermi surface. 
The pairing susceptibility is computed only for the most 
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FIG. 5. Projection of momenta on the Fermi surface; dis- 
cretization and labelling of angle variables. 

interesting case of vanishing total momentum q = 0. The 
corresponding vertex part is approximated by 

i?^(0;k,-k) « i?^(0;kj.,-kj.) (56) 

The projection approximation for the vertex parts is 
again exact for A = Aq and asymptotically exact for 
A — > 0. It can also be improved for intermediate A by 
projecting momenta either on the Fermi surface or on the 
van Hove surface, whichever is closer. 

The behavior of the compressibility xc = dn/dfi = 
Xc'(O) and the homogeneous spin susceptibility xs — 
Xs(0) cannot be obtained directly from the flow equa- 
tions for Xc,s(l)- "^^^ problem is that Xc,si^) vanishes 
for all A > (at zero temperature), since the infrared 
cutoff blocks particle-hole excitations with an infinitesi- 
mal momentum transfer. We will therefore compute xc 
and xs from the effective quasi-particle interaction 

/kTkt - ^k.< r^kpa, k'^a'; k^a, k'^a') (57) 

Note that the forward scattering limit (zero momentum 
and energy transfer) of the two-particle vertex is unique 
for A > and converges to the quasi-particle interaction 
for A ^ and kp k'p |2^]. The wave function renor- 
malization factor is one in our calculation because 
we have neglected self-energy contributions. Following 
the usual Fermi liquid arguments p8| one obtains the 
compressibility and the homogeneous spin susceptibility 
as 

Xcs = 2 J ^ X^f 6(6^ - ^l) (58) 
c s ' 

where X^^^ is the solution of the inhomogeneous linear 
integral equation 

X^f + 2/1^ f^;i, X-f Sie., - M) - 1 (59) 
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h (/k^k'^ ± fZ^) ■ In a non-interacting 



with fF'^ , 

system one would obtain X^'^ = 1. The quantity 



dfj, 



1 



Vk^ 



(60) 



with the velocity Vk — VkCk describes the linear response 
of the Fermi surface in (a shift along the normal vec- 
tor) for a small shift of the chemical potential. 

For a concrete numerical solution of the flow equa- 
tions the angular dependence of the vertex function is 
discretized, with a finer mesh in the vicinity of the saddle 
points of Ek at (tt, 0) and (0, tt), as shown in Fig. 5. One is 
thus left with a finite (though very large) number of flow- 
ing coupling constants. If not stated otherwise, we have 
used 16 points on the Fermi surface to discretize the ver- 
tex function, corresponding to 672 independent (i.e. not 
symmetry-related) coupling constants, if the fermi level 
lies at the van Hove energy e^HS, and 846 couplings oth- 
erwise. With the addition of 16 points on the van Hove 
surface one has to deal with 4728 flowing couplings. 



C. Results 

We have computed the flow of the vertex function and 
the susceptibilities for several choices of the bare interac- 
tion U > 0, the next-nearest neighbor hopping amplitude 
t' < and the chemical potential fj,, where t' and fi have 
been chosen such that the Fermi surface is on or close 
to the van Hove points of ek, and the particle density is 
close to half-filling. 

In all cases the vertex function develops a strong mo- 
mentum dependence for small A with divergencies for 
several momenta at some critical scale Ac > 0, which 
vanishes exponentially ?7 — > 0. To see which physical in- 
stability is associated with the diverging vertex function 
we have computed the following susceptibilities: 

i) commensurate antiferromagnetic spin susceptibil- 
ity Xs(7r,7r), 

ii) incommensurate spin susceptibility xs(q) with q = 
(7r-(5,7r) andq= (1 - (5)(7r, tt) @, 

iii) commensurate charge susceptibility xciT^TT^) j 

iv) singlet pair susceptibilities with form factors |^ 

1 (s-wave) 
I -^(cosfca; + cosfcj^) (extended s-wave) 
' -^{cos kx — cos ky) (d- wave 1^2.2 _j^2) 
(d-wave dxy)- 

(61) 



sin kx sin ky 



Some of these susceptibilities diverge together with the 
vertex function at the scale Ac. Depending on the choice 
of [/, t' and the strongest divergence is found for the 



commensurate or incommensurate spin susceptibility or 
for the pair susceptibility with dx^-y-^ symmetry. 

We will now present explicit results for the flow of 
the two-particle vertex and susceptibilities for a coupling 
strength U — t, which is much smaller than the band- 
width W — 8t and therefore safely in the weak coupling 
regime. All energy scales will be plotted in units of t. 
To exhibit the interaction-induced renormalizations of 
the susceptibilities, we plot the flow of the ratio X^/Xo' 
where Xo is the susceptibility of the non-interacting sys- 
tem at scale A, as obtained from the flow equations for 
U = 0. We show examples for the flow of Xo in the Ap- 
pendix. Note that the non-interacting susceptibilities Xo 
are all finite for A > 0, such that a divergence of x^ at a 
finite scale Ac implies a diverging ratio x^ /Xo and vice 
versa. 

In Fig. 6 we show the flow for t' = and fi = —0.005, 
corresponding to a density n — 0.995, i.e. almost at half- 
filling. Here and in the following we plot the singlet part 
of the vertex function for a selected choice of momenta on 
the Fermi surface, including those momenta for which the 
vertex function renormalizes most strongly. The singlet 
vertex function has its largest values for umklapp scat- 
tering along the diagonal of the Brillouin zone, but also 
forward and Cooper scattering of particles on opposite 
sides of the almost square Fermi surface are strongly en- 
hanced. Scattering amplitudes for momenta near the van 
Hove points diverge a bit more slowly. The triplet part 
of the vertex function is renormalized mostly for forward 
and Cooper scattering, but generally more weakly than 
the singlet part. The spin susceptibility with an antifer- 
romagnetic wave vector clearly dominates over pairing 
susceptibilities in this case. The incommensurate spin 
susceptibilities are indistinguishable from the commen- 
surate one in Fig. 6 because the incommensurability pa- 
rameter S is almost zero so close to half-filling (see Ref. 
|p9| ) . Note also that the susceptibility ratios for isotropic 
and extended s-wave pairing are equal here, and almost 
coincide with the charge density susceptibility ratio. The 
non-interacting susceptibility for extended s-wave pairing 
(and thus x) is however much smaller than the other two 
(see Appendix). 

Decreasing the density (away from half-filling) one en- 
ters a regime where pairing correlations with dj.2_j,2 sym- 
metry dominate at sufhciently low energy scales. This is 
seen in Fig. 7, where we show the flow for t' — and 
/i = —0.02, corresponding to n = 0.984. Note that for 
small U the transition from antiferromagnetism to super- 
conductivity occurs already at a critical density ric quite 
close to half-filling. For increasing U the deviation of ric 
from half-filling increases. The flow in Fig. 7 exhibits a 
threshold at A = 2|/i| below which the amplitudes for var- 
ious scattering processes, especially umklapp scattering, 
renormalize only very slowly. The flow of the antiferro- 
magnetic spin susceptibility is cut off at the same scale. 
The infinite slope singularity in some of the flow curves 
at scale A = |^| is due to the van Hove singularity being 
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FIG. 6. (a) The flow of the singlet vertex function as 
a function of A for several choices of the momenta Rfi, k_F2 
and \i'pi, which are labelled according to the numbers in Fig. 
5. The model parameters are U = t and t' = 0, the chemical 
potential fi — —0.005; (b) the flow of the ratio of interact- 
ing and non- interacting susceptibilities, IXo for the same 
system. 

crossed at that scale. The pairing susceptibility with 
rf3;2„j,2-symmetry is obviously dominant here (note the 
logarithmic scale). Following the flow of the ver- 
tex function and susceptibilities one can see that 
the o?2,2_j,2-pairing correlations develop in the presence 
of pronounced but short-range antiferromagnetic spin- 
correlations, in agreement with earlier ideas on d-wave 
superconductivity 

In Fig. 8 we show the U) phase diagram for t' — 
obtained by identifying the dominant instability from the 
flow for many different values of /i and U. Note that 
fi — corresponds to half-filling. The regime with a 
leading commensurate antiferromagnetic spin density in- 
stability is separated from the d-wave pairing regime by 
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FIG. 7. Same as Fig. 6 ioi U = t, t' = and ^^ = -0.02. 

a thin region where incommensurate spin density fluc- 
tuations with q = (tt, TT — (5) dominate. Other incom- 
mensurate structures may be more favorable than the 
ones considered here. For [/ — *■ at flxed density n < 1 
the superconducting instability always dominates, be- 
cause the bare particle-hole bubbles are flnite away from 
half-filling, while the Cooper channel always diverges at 
least logarithmically. The way the critical energy scale 
Ac varies as the system is doped away from half-flUing 
can be seen in Fig. 9 for an interaction strength U ~ t. 

The different symbols show which instability is lead- 
ing at Ac. The two straight lines represent the linear 
functions Ac = and Ac = respectively. As al- 

ready observed by Zanchi and Schulz ||l^], the supercon- 
ducting instability is leading if Ac < This may be 
related to the fact that only pair fluctuations receive a 
singular enhancement at A = |/Lt| (see Fig. 7), while spin 
fluctuations don't. A commensurate spin density wave 
instability cannot be favorable for Ac < 2|/x|, since their 
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spin-density wave regime from the superconducting regime, 
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surate spin-density regions. 
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FIG. 9. The critical energy scale Ac as a function of the 

chemical potential fi for U = t and t' = 0. The different 
symbols indicate, whether the leading instability is a com- 
mensurate or incommensurate spin-density wave or d-wave 
pairing instability; the straight lines represent the functions 
Ac = and Ac = 2|/i|, respectively. 

flow is cut off at A = 2|/i| (see Fig. 7 once again). Hence, 
the incommensurate spin density wave is the leading for 
/i| < Ac < 2|/i|. The sharp peak in A,; at /i = (half- 
filling) is due to the van Hove singularity. For larger 
deviations from half-filling, the critical energy scale Ac 
vanishes rapidly. Note, however, that for larger values of 
U the regime with a sizable scale Ac extends to larger 
values of n, i.e. to larger doping. 

In Fig. 10 (a) and (b) we show results for the compress- 
ibility and the homogeneous spin susceptibility, respec- 
tively, for two choices of the chemical potential, atU = t 
and t' = 0. We recall that these quantities have been ob- 
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FIG. 10. The flow of (a) the compressibility and (b) the 
homogeneous spin susceptibility as a function of A for vari- 
ous choices of /I at t/ = 1 and = 0; and x"'" are the 
corresponding non-interacting quantities. 

tained from the forward scattering vertex by using Fermi 
liquid relations as discussed above. The non-interacting 
compressibility kP and spin susceptibility x*'° in the plot- 
ted ratios are defined without infrared cutoff. Hence the 
flow in Fig. 10 is entirely due to the flow of the Landau 
function, starting at the simple RPA result for the Hub- 
bard model at A = Aq. Close to half-filling, where a spin- 
density wave instability is leading, the compressibility is 
suppressed at low energy scales, as expected for a system 
with a charge gap at or near the chemical potential. The 
homogeneous spin susceptibility remains finite near the 
spin-density wave instability. By contrast, further away 
from half-filling in the regime where the d-wave pairing 
instability is leading, the compressibility diverges while 
the homogeneous spin susceptibility is suppressed. A 
suppressed spin susceptibility is expected as a precursor 
of the spin gap opening in any spin singlet superconduc- 
tor. Very close to the instability the spin susceptibility 
flows through zero to negative values, which implies that 
our 1-loop calculation breaks down in this strong cou- 
pling regime. A diverging compressibility would indicate 
a tendency towards phase separation, but the increase of 
K sets in quite abruptly only very close to the instabil- 
ity, where the renormalized couplings are already so large 
that the 1-loop results are not reliable any more. In any 
case the large charge fluctuations indicated by k — > oo 
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FIG. 12. Same as Fig. 6 for [/ = t' = -0.05 and = 4t' . 



would only be a consequence of the pairing instability of 
the Hubbard model, not a driving mechanism, since the 
pairing correlations appear already at a higher energy 
scale. 

Results for the flow of the vertex function and sus- 
ceptibilities for t' < and fi — e^n — 4t' are shown in 
Figs. 11 and 12, with t' = -0.01 and t' = -0.05, re- 
spectively. The corresponding Fermi surfaces touch the 
saddle points at (7r,0) and (0,7r). In the first case the 
density is n = 0.992 and in the second one n — 0.959. 
For the bare interaction we have chosen U = t as before. 
The major difference with respect to the perfect nesting 
case t' = is that now the umklapp processes near the 
diagonal of the Brillouin zone are much less enhanced 
at low energy scales, such that scattering processes with 
momenta near the van Hove points (tt, 0) and (0, tt) be- 
come most prominent. In this situation the simple scal- 
ing approaches which concentrated exclusively on the van 
Hove points ||l^,|l3| provide already a useful qualitative 
picture of the important effective interactions and their 



renormalization. Antiferromagnetic correlations are now 
mostly driven by umklapp processes from (tt, 0) to (0, tt) 
and vice versa which, due to the equivalence of the points 
(tt, 0) = (-TT, 0) and (0,7r) = (0, -tt) in the Brillouin 
zone, can also be viewed as Cooper processes. Indeed 
these processes are also responsible for d-wave pairing 
correlations. 

For the parameters chosen in Fig. 11 antiferromagnetic 
correlations dominate over pairing. The incommensurate 
susceptibility with q = (tt, tt — (5) is a bit larger than the 
other incommensurate candidate and the commensurate 
antiferromagnetic susceptibility. Note that there may be 
other still larger incommensurate susceptibilities among 
those not computed here. We have merely investigated 
two (frequently discussed) incommensurate spin suscep- 
tibilities out of a variety of infinitely many candidates. 
Moving further away from half-filling, as in Fig. 12, one 
finds again dominant pairing susceptibilities, with (i^2_j,2 
symmetry in each case. 

The phase diagram in the (<', J7)-plane with /i = 4<' < 
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FIG. 14. The critical energy scale Ac as a function of t' for 
= 4t' < and U — t; the short dotted hne separates the 
spin density regime from the d-wave pairing regime. 

is shown in Fig. 13. Note that the chemical potential is 
always situated at the van Hove singularity here and the 
density decreases away from half-filling with increasing 
\t'\. Since we have no good guess for the optimal density 
dependence of the incommensurability vector for t' ^ 
we have not distinguished different spin density waves in 
Fig. 13. The behavior of as a function of < with 
/i — At' and U — t is shown in Fig. 14. 

The decrease of Ac with increasing (and thus in- 
creasing doping) is slower here than in Fig. 9, since the 
Fermi level remains on the van Hove singularity such that 
only the importance of nesting is weakened under doping. 

All the numerical results discussed above have been ob- 
tained by projecting momentum variables of the vertex 
function on 16 points on the Fermi surface as shown in 
Fig. 5. To see how much these results may be modified in 
a more refined projection scheme, we have computed the 
flow for some typical model parameters with a projection 




FIG. 15. The critical energy scale Ac as a function of the 
number of discretization points A'^o on the Fermi surface, for 
U = t,t' = -O.OU and fi ^ 4t' . 

on 32 points on the Fermi surface, and also with a projec- 
tion on 16 Fermi surface points and 16 additional points 
on the van Hove surface. It turned out that these refine- 
ments, which increase the computational effort consider- 
ablely, lead only to a moderate reduction of the critical 
energy scale, without changing the qualitative behavior 
of the vertex function and susceptibilities. In Fig. 15 we 
show the dependence of the critical scale Ac as a function 
of the inverse number of discretization points A^o on the 
Fermi surface for A'o — 4, 8, 16, 32 and a fixed choice of 
model parameters. We see that the critical energy scale 
obtained from a discretization with 16 points has already 
the right order of magnitude. 



IV. CONCLUSION 

In summary, we have shown that the renormalization 
group method developed by Salmhofer with our ex- 
tension for the computation of susceptibilities can be 
used as a systematic tool for detecting instabilities in 
a weakly interacting Fermi system with several coupled 
infrared singularities. Such an RG analysis is completely 
unbiased. The selection of retained Feynman diagrams 
is dictated by the weak coupling expansion and can be 
systematically improved by including higher orders in a 
loop expansion. 

Evaluating the fiow equations on 1-loop level for the 2D 
Hubbard model we have found antiferromagnetic insta- 
bilities close to half-filling and dominant superconducting 
instabilities with dr^2_y2 symmetry at smaller densities 
(still near half-filling). Incommensurate spin structures 
can be favorable in the antiferromagnetic regime near 
half-filling. 

The critical energy scale Ac where vertex functions and 
susceptibilities diverge vanishes exponentially as C/ — > 0, 
but becomes sizable already for relatively weak coupling 
strengths (compared to the band width), even in the su- 
perconducting regime. The appearance of strong pair- 
ing correlations with d^2_y2 symmetry in the 2D Hub- 
bard model at physically interesting energy scales is thus 
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well established at weak coupling. The flow of the vertex 
function and susceptibilities clearly shows that the pair- 
ing instability is driven by short-range antiferromagnetic 
correlations in the system. This supports earlier ideas 
and numerical results (for finite systems) suggesting d- 
wave superconductivity driven by antiferromagnetic cor- 
relations in the Hubbard model Note that Ac must 
not be interpreted as a transition temperature for antifer- 
romagnetism or superconductivity, but rather as an en- 
ergy scale where bound particle-particle or particle-hole 
pairs are formed. A Kosterlitz-Thouless transition to a 
superconducting state may occur at a lower energy scale 
while antiferromagnetic order is of course possible only 
in the ground state of a two dimensional system with 
spin-rotation invariance. 

We finally outline some interesting extensions of the 
present work for the future: 

i) Non-local interactions: Non-local interactions may 
play an important role even though they are usually much 
smaller than the local (Hubbard) interaction. They af- 
fect the RG flow via a different initial condition for the 
vertex function and can thus be taken into account very 
easily. 

ii) Fermi surface instabilities: The Fermi surface is gen- 
erally deformed by interactions. Computing a suscepti- 
bility for Fermi surface deformations from the RG flow 
one finds that deformations breaking the discrete square 
lattice symmetry may occur Q. 

iii) Self-energy effects: It will be interesting to compute 
self-energy contributions and see how they affect the in- 
stabilities. The numerical effort for this is small on 1-loop 
level, but also a 2-loop calculation seems feasible. Kishine 
and Yonemitsu [^l| have recently computed the renor- 
malization of the quasi-particle weight on 2-loop level for 
two flat Fermi surface pieces, but the feedback of self- 
energy effects on instabilities has not yet been treated. 
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APPENDIX A: NON-INTERACTING 
SUSCEPTIBILITIES 



Here we show results for the flow of the non- interacting 
susceptibilities Xo for the choices of t' and fi correspond- 
ing to those in Fig. 6, 7, 11 and 12. The reader may thus 
estimate the absolute scale of by multiplying Xo '^ith 
the results for the ratios x^/Xo '^^ -'^^I- 
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FIG. 16. Free susceptibilities for (a) t' = and 
/i = -0.005f, (b) = and /i = -0.02f, (c) t' = -O.Olt 
and fi = 4:t' and (d) for t' — —0.05t and jj, — 4t' , correspond- 
ing to the examples in Fig. 6, 7, 11 and 12, respectively. 
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The incommensurate spin-density susceptibilities and 
the charge-density susceptibiUty lie too close together to 
be always individually seen. The extended s-wave pair- 
ing susceptibility is of the order of 10~^ and therefore out 
of scale. 
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